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l. INTRODUCTION 


The purpose of this paper is to investigate and recommend suitable 
types of measurements on the lunar surface which will result in the 
establishment of a selenodetic control network. In this paper, we will 
define a selenocentric celestial coordinate system, in which we may 
make measurements, and in turn determine an astronomic position on 
the lunar surface. Combining this astronomic position with a knowledge 
of the deflection of the vertical at this position, we could then determine 
the coordinate of this position in a selenographic system. 

Therefore this paper will be divided into two primary areas which 
are: 

(1) Establishment of a selenocentric celestial coordinate system 

and measurements in that system from the lunar surface. 

(2) Determination of the deflections of the vertical on the lunar 


surface. 








2. ESTABLISHMENT OF AND MEASUREMENTS IN THE LUNAR 
ASTRONOMICAL COORDINATE SYSTEMS 


On the earth most of the geodetic astronomic calculations are made 
using a geoequatorial system where the star coordinates are right ascen- 
sion and declination. There are a great number of catalogs, ‘maps and 
almanacs giving the coordinates of the stars and other celestial bodies 
in this system. If the coordinates in this system could be easily trans- 
formed to a lunar equatorial system, all the methods of geodetic astronomy 
could be employed on the moon for selenodetic purposes. In this section 
we will show such a transformation complete with error analysis, while 
also discussing such related factors as time, location of lunar pole stars 
and observations in the horizon system. 

We will also show some methods which may be used on the lunar 
surface to determine astronomic position. These methods will be dis- 
cussed in relation to the instrument requirements, and accuracy of the 
system as a whole. Included among these methods is an inertial system, 
which is presently installed on the lunar module, and which in itself can 
establish a position on the lunar surface. 

The methods and instruments presented are not assumed to be a 
complete listing, but are based on practical factors which may obviously 
limit other methods because of physical size or complexity of equipment, 
time required to achieve an astronomic position, or excessive burden 
being placed on the astronauts, either during training or during the 
operation itself. In this light we would most like to find a method which 
is as conservative of the astronauts time as possible (in both training and 
employment), achieves satisfactory results, and does not involve undue 
calculations or recording procedures. In conjunction, we wish the 
instruments involved in this method to be simple, portable, and of 


e 
_— 


course accurate. 





2.1 Lunar Astronomy 


Kolaczek [14] defines a selenoequatorial coordinate system (lunar 
right ascension system) where the axis of rotation of the moon defines 
the celestial pole, which is the basic direction, and the celestial 
equator, which is a plane perpendicular to the rotation axis, containing 
the center. The latter is defined as the basic plane. This system is 
affected by precession and nutation (physical libration) of the lunar axis 


of rotation. The selenoequatorial system is analogous to the geoequatorial 


system on earth, but the selenoequatorial coordinates of the stars change 
their values faster than the geoequatorial coordinates due to the lunar 
precession being 1,360 times faster than the earth's. The selenoequatorial 


coordinates are: 


@ - Lunar right ascension measured on the lunar celestial 
equator from its ascending node on the ecliptic. It is 
measured to the east from 0° to 360°. 

d - Lunar declination which is the angular distance from the 
lunar celestial equator. It is measured from 0° to 90°, 
positive in the northern celestial hemisphere, and negative 
in the southern. 

Now we will take a coordinate system with known coordinates and 
transform these coordinates into the selenoequatorial coordinate system. 
The translation of any coordinate system from one point in space to 
another changes the values of the spherical coordinates of fixed points 
on the celestial sphere. This change is caused by the parallax or the 
translation of the origin of a coordinate system and the aberration caused 
by the different motion of this newly translated system. The translation 
of mean earth oriented coordinates into mean lunar oriented coordinates 
does not require consideration of parallax if we consider only the coordi- 
nates of stars, aS we may omit the parallax of the earth-moon distance, 
due to the great distance of the stars, as compared to the earth-moon 
distance. The transformation of the apparent coordinates however requires 


consideration of the influence of the moon's aberrations. 
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2.1.1 Transformation of Mean Geoequatorial into Mean Selenoequatorial 


Coordinates 

The transformation of mean eeoequatorial into mean selenoequatorial 
coordinates can be performed in two steps []4]. First the mean geoequa- 
torial coordinates &@, dare transformed into the mean ecliptic coordinates 
X, Bby the well known formulas: | 

sinB = cos€ sin6 - sin€ cos sina 

cosB cosA = cos cos@ 

cos sind = sin€ sin6 + cos€ cos6 sin@ 
where € the obliquity of the ecliptic is 

€ = 23°.452294 - 0 .0130125 T - 0°.00000164 T* + 0°.000000503 1 
where T denotes the time measured in Julian centuries of 36,525 ephemeris 
days from the epoch 1900 January 0.5ET = J.D. 2415020.0. 

Next the transformation formulas of mean ecliptic coordinates into mean 


selenoequatorial coordinates. 


sindz = sinB cosI + cos8 sinI sin( A - §) 
cos dx ame = sin8 sinI - cosB cosI sin (A - ©) Ziel 
cos dz cosa? = - cosB cos(A - &) 


where I is inclination of the lunar equator to the ecliptic 
I= 1° 32' o1" 
{2 The mean longitude of the ascending node of the lunar orbit on 
the ecliptic is 
Q= 259°.18305556 - 1934°.14200833 T + 0°.002078T* + 0.000002 T® 


where T is defined as above. 


2.1.2 Transformation of Mean Selenoequatorial Into Apparent Selenoequatorial 


~ 


Coordinates 
The transformation of mean into apparent coordinates is given in 


a form [14] similar to that used for the earth's equatorial coordinates 


¢2 QQ 7 
aap a ne eee io On be Olney 1) d, + Ta a 
Hlecappe ft eC ete Tete ois, ch Cyeey + Dice Ta 114 
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Similarly to Bessel's equations on the earth, Ai and B; are the effects 


of precession and nutation to date from the determination of the mean 


coordinates. The Ci and Di terms are the effect of lunar annual abberation, 


ane Tt. 1 (4/2 ) terms are the corrections for proper motion. The 


ai,b1,01,d1,23,bi,cl,di coefficients are the lunar Besselian star constants. 


Specifically: 

Ai = -(Pet + @)sinI 

1S) eee 

Ci = - K sin(Iap - 2+ 180°) cosI 

Di = - K cos(Lap - + 180) 

where: 

De Daily motion of the ascending node 
of the moon's orbit on the ecliptic 
(= - 0°. 0529539222) 

t Number of ephermeris days from the 
epoch t. of the mean coordinates @ 2 da 

I Inclination of the moon's equator to the 
ecliptic (=1° 32’ 01") 

PO Physical librations in inclination and in 

= ascending node respectively (see Appendix #3) 

K,k Constant of the earth's and of the lunar 
aberration respectively . 
K = k(1 - Vi cos(Le - 1:)/ Ve) 
(k = 20'. 4958) 

Ls,li; Lap Longitude of the sun, the moon and the 


moon's apex in the motion around the 


sun respectively 


lae= be 0) ce Aa 





7 = 
re 





) Vi sin (L, - h ) 
A= i : 
areas (2 ik Lt | 


Vigra - Velocity of the earth's and the moon's 


orbital motion, respectively 


l| 


29.8 km/sec 
1.02 km/sec) 


(Ve 
Var 


The lunar Besselian star constants are: 


I 


QQ 
al (cotI - sin @, tandn) 
al — COS Ol. 
QQ 
bi = cos @ , tan dy 
bz =. sin & i] 
Q2 
Ci =- cos@, sec dy 
' Pee 
Gq =sin@, sin dz + tan I cos ds 
_ a 
di =sin@, sec dz 
i Q 
di =cos @, sin dy 
The proper motion will take on similar values to those computed 


on the earth: 


Ta t 


Ha, Wa = The components of proper motion of the star 


\ 


2.1.3 Lunar Horizon System 


Measurements on the lunar surface will most likely be made in the 
horizon system. In the moon's horizon system, altitude ain and azimuth Aan 
are defined in the same way as that of the earth's horizon system. The 
transformation between the lunar horizon system and the lunar selenoequatorial 


system similar to that given by Mueller [16] is as follows: 


QQ 
cos d sin ( LHA,-.A-@ ) =- sin Z: sin Aan 


2 


sin d = cos Z1 sin®+ sin 27 cos A cos? 222 


QQ 
cos d cos (LHA, -X-@ ) = cos Zi cos@- sin Zi cos Ai sin® 
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where LHA - The lunar hour angle of the node measured from the 
a 


reference median 
© - The lunar astronomical latitude 
X - The lunar astronomical longitude 


ZZ. = 90"— ain 


The above equations can be derived with the use of Figure 1 below: 


NCP 


90-d 
{2 


h = LHA, -r’ -a@ 





Z1 =90-arp star 


Figure 1. - The Astronomical Triangle 


Zae amar lime 


As we have defined the lunar right ascension system in terms of 
ascending node, it would only be appropriate, that we define the lunar time 
system in terms of the draconic month. The draconic month is defined by 
two successive passages of the moon through the same node of its orbit. 
The draconic month is equal to 27.212220 mean solar days. 


The expression of selenographic longitude in terms of time would be 


as follows: 
imeanionth  —  2%.2'2220 mean solar days = 360° 
1 day = 9070740 mean solar days = 12° 
1] hour = .03779475 mean solar days = _ 30! 
1 minute = .0006299125 mean solar days = 380" 
1 second = .0000104985 mean solar days = Ome 








Figure 2 below shows the relationship in draconic times. 









eIpIqowl oP 


local npean. a 
astronamic meridian 


LMDT Lunar mean draconic time 

LADT Lunar apparent draconic time 
RMDT Reference mean draconic time 
RADT Reference apparent draconic time 
eq.N Equation of the node 


Figure 2. - Draconic Time System 


From Figure 2 it is evident that: 
LMDT + eq.N = LADT 
and LMDT - RMDT = A 


where A - The reduced astronomic longitude with the positive 
sign to the east 
and eq.N - The lunar equation of the node which is due to nutation 
| (physical libration in longitude T and in node 0) (see 


Appendix C). 


_— 





From the actual rotation angles of the moon, we find that eq.N =T-o. 


We will assume in this section that the meridian is not affected by 


polar motion and that the rotation of the moon is free from short periodic 
irregularities. 


The hour angle of the node is equal to the lunar right ascension of 
an object at the instant of its upper transit over the meridian from 
which the hour angle is reckoned. Therefore, draconic time can be 


determined by observing transits of stars. 


The lunar time intervals are almost equal to the terrestrial time 
intervals of the same name (91% of the terrestrial values), which is 
still advantageous from the psychological point of view because a change 
to such a time system would be almost unnoticeable to man. The 


meridian to use for a change of dates could be the 180° meridian, 


2.3 Lunar Pole Stars 


As stated by Gurevich [5] the lunar poles describe circles in the 
celestial sphere with an angular diameter of 3°. This circle is centered 
on the ecliptic pole, with a period of {18.6 years. Therefore, first 
certain stars, then others, will be most suited for use as pole stars. 
Tables of altitudes and azimuths must be compiled for the stars closest 
to the poles at a particular time. The mean north pole 
of the moon described in the earth right ascension system would be 
approximately, @ = 18), 6 = 66°33'. We see from the FK4 catalogue 


that the following stars from Table 1 are near the north pole. 





s 








Table 1. Lunar North Pole Stars 


Catalog Number Magnitude a 5 
659 5.21 le ee" 68° 09! 
701 6.00 18° 36" 65° 27! 
664 4.87 17° 37" 68° 46! 
G85 | 5.03 I 12 64° 23! 
and approximately @ = Ga , O=F- 66° 33' for the south pole in Table 2 


from the Boss catalogue. 


Table 2. Lunar South Pole Stars 


Catalog Number Magnitude O 6 
6927 5.3 5 32" - 64° 15" 
7246 4.5 5" 44" ~ 65° 45! 
7384 5.2 5° 49" ~ 66°55! 
7813 5.8 6° 06" ~ 66 02! 
7946 4.9 6° 11" - 65° 35! 


With a knowledge of the relationship of the position of the lunar pole to 
those stars the astronauts can make a rough north alignment with a probable 
accuracy of less than 1 minute while using only a gun sight type sighting 


device. 
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2.4 Measurements on the Lunar Surface 


Presented here are four basic types of measurements which can be 
employed on the lunar surface to determine astronomic position. 

One method is the measurement of altitudes or zenith angles of at least 
two stars separated in azimuth, combined with the known celestial coordinates 
of a star and a knowledge of time. This is a method which can be employed 
with an ordinary sextant or theodolite. 

A second method is by determination of a relative bearing and altitude of 
a star combined with information from inertial and gravitational sensors, which 
is the method used by the Alignment Optical Telescope (AOT) on the lunar 
module. 

A third method is by measurement of relative azimuth angles between at 
least three known stars. This technique can be used with an ordinary theodolite, 
but also lends itself well to the use of an automatic instrument. 

A fourth method is by star photography, which can use either the first 
or third, or other methods to determine position. Here we can use a photo 
theodolite, or other metric camera. 

In the following section these four methods will be discussed in relation 
to their ease of employment on the lunar surface, instrumental adaptability 


and accuracy, and other pertinent facts. 
2.4.1 Method 1: Measurements of Altitudes 


The measurement of altitudes or zenith angles in order to determine 
astronomical positions can be accomplished by the use of a sextant or 
theodolite. The obvious advantage of this method lies in the fact that both 
sextants and theodolites are comparatively light equipment, which could 
feasibly be included on a lunar mission. The sextant is a familiar 
instrument to the astronauts, as they already have been trained in its use. 
Although, the training required to allow the astronauts to effectively use a 


theodolite could be accomplished in a short period of time. The big 
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advantage in using this equipment is that astronomical positions can be taken 
at several points. 

As stated in the Geonautics study [7], it is not likely though that a 
theodolite would be practicable for the initial missions. The time and effort 
required to set up and measure is far out of proportion to the many other 
tasks required of the astronauts. It would also be necessary to adapt the 
fine pointing and adjusting screws for use with the pressure gloves the 
operator must wear. In addition, precise sightings through the helmet 
face plate would be most difficult. 

A hand held sextant may be used to determine a position on the lunar 
surface. The altitude of a star is read directly from the scale. With the 
_ aid of sub-divisions and a vernier, altitudes can be read with a first-class 
instrument to one-tenth of a minute of arc. 

The problem which arises with the use of a simple sextant is in the 
determination of the lunar horizon. In order to eliminate this problem, 
we can consider the use of an artificial horizon or a bubble sextant, so 
that we then need only define the local vertical instead of the horizon. 

The artificial horizon is a rectangular shallow basin of mercury. The 
observer places himself so that he can see the image of the body, whose 
altitude is to be measured, reflected in the mercury. He points the 
telescope toward the artificial horizon so that he uses a second image of 
the heavenly body in place of the horizon. By moving the index arm, the 
two images are made to coincide. The circle reading is twice the altitude 
which thereby reduces the precision of the instrument to 12 seconds of are. 

A bubble sextant (bubble octant) [17] contains its own artificial horizon. 
In measuring altitudes, the image of the celestial body is brought into 
coincidence with the image of a circular level bubble which forms the 
artificial horizon. The altitude is then read on a micrometer drum. ‘The 
fine drum is divided into 2 minute intervals which can be interpolated to 
the nearest 1 minute. In a-bubble sextant such as the Bendix ''Pioneer", 
provision is made to record a series of altitudes automatically. An 


average reading can be determined from these readings. 
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These instruments are not capable of first order geodetic accuracy. 
The positional error can be a maximum of 50 meters from the 
sextant with the artificial horizon, and 600 meters for the bubble sextant. 
This, of course, does not take into account, sextant internal errors, 
bubble errors, and errors in the star tables. In addition, the recording 
of each measurement using the common sextant will burden the astronauts. 
With both instruments each star must be identified by the astronauts, 
leading to further possibilities of error. Time recording must also be 
done manually. 

Using this method (see Figure 1) we can find azimuth from the 


formula as given by Mueller [16] 


- cos d sinh 


sin A = : 
sin Z 


where h = LHA g- A -— g& 


We can find the latitude from the basic formula as given by Mueller [16] 
cosZ = sind sin@+cosd cosh cos® 
and we can determine longitude from the formula as given by Mueller [16] 


cos Z - sind sin® 


cos h = 
cos d cos® 


2.4.2 Method 2: Measurement with the Alignment Optical Telescope 


A relative bearing and altitude of a star can be determined by the 
Alignment Optical Telescope (AOT) which is currently installed in the Lunar 
Module (LM). An astronomic position on the lunar surface can be established 
by the use of this instrument in conjunction with the Inertial Measurement 
Unit (IMU) and the Primary Guidance Navigation System (PGNCS) accel- 
erometers. The great advantage of this method lies in the fact that no 
additional equipment or training of the astronauts is necessary to obtain an 
astronomical postion. An apparent disadvantage is the fact that only one 


astronomical position can be determined during each landing due to the 
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fact that the equipment is fixed in the lunar module which in turn is fixed 
on the lunar surface. 

The AOT as described in the Grumman manual [4] is a unity power 
periscope type device with a 60° conical field of view. Its primary pur- 
pose is to provide a method of establishing an inertial reference with 
respect to celestial bodies for inflight and lunar surface operations. It 
is operated manually by the astronauts and is mounted on the PGNCS 
navigation base to establish a mechanical alignment and common reference 
with the IMU. The AOT has a moveable shait axis which is parallel to 
the LM ''X" axis and a Line of Sight (LOS) axis which is the center of 
field of view and fixed at 45° from the LM + "X" axis. 

The AOT LOS is fixed in elevation and moveable in azimuth to six 
detent positions. These detent positions are located at 60° and are 
selected manually by turning a detent selector knob on the AOT. Three 
of the positions can be used operationally for star sightings and are called 
forward, left, and right (see Figure 3). The forward, or zero detent 
position places the LOS in the x,z plane looking forward and up. The 
rieht, or +60° position places the LOS 60° to the right of the x,z plane and 
the left or -G0° position places the LOS 60° to the left of the X,Z plane. 
Each of these positions maintain the LOS at 45° from the LM + X axis. 
The remaining three detent positions are storage positions with the 
position 180° from the forward position being the normal storage position. 
The reticle control makes it possible for manual rotation of the reticle 
for use in lunar surfece alignments. The AOT reticle pattern consists of 
| cross hairs and superimposed archimedes spiral (see Figure 4). 

The vertical cross hair is an orientation line designated the 'Y" line 
and is parallel to the LM X axis when the reticle is at the 0 reference 
position. The horizontal cross hair is an auxiliary line designated the 'X" 
line and is perpendicular to the orientation line. A one turn spiral is 


Superimposed from the center of the field of view to the top of the vertical 
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FIGURE 3. - AOT Line of Sight Positions 


15 





SNIT NOILYINS WO 
OL WRYON 
SMITA KAY 


xX¢ 


\ 








uLayeg eTONeY LOV - 


WuldS 


3INSY3 438 
OW3Z 


° ? GdNolLt 


009 MBIA 490 OlAala 


2NEV NOLLYLN3ATHO 


16 





cross hair. The reticle is edge illuminated. 

A counter located on the left side of the AOT, displays an angular 
readout of reticle rotation. The counter reads in degrees to within +0.02 
degrees or +72 seconds. Interpolation is possible to within +0.01 degrees. 

On the lunar surface, the selected star can be sighted in either the 
AOT left,right or forward positions. The astronaut using the manual drive 
knob, adjusts the reticle to superimpose the orientation line on the target 
star. This reticle angle displayed on the AOT counter is then manually 
inserted into the Lunar Module Guidance Computer (LGC) via the Display 
Console (DSKY). This provides the computer with the star orientation 
angle (shaft angle). The astronaut then continues to rotate the reticle until 
a point on the spiral is superimposed on the target star. This second 
angular readout is then entered into the LGC via the DSKY. The AOT 
detent position and the star code number are also inserted into the LGC 
via the DSKY. The LGC can now calcuiate the angular displacement of the 
star from the center of the field of view by computing the difference between 
the two counter readings. Due to the characteristics of the reticle spiral, 
this Aangle is proportional to the distance of the star from the center of 
the field of view. Using this Aangle and the proportionality equation. the 
computer is then able to calculate the trunnion angle. During the lunar 
surface phase of the lunar landing mission, the Mission Control Center 
and the LGC use the unit gravity vector as measured by the LM PGNCS 
accelerometers to align the LM platform and to estimate the latitude and 


longitude of the LM on the lunar surface. 


The AOT star alignment establishes an inertial reference by employing 


the following formulas: 


> 5 = > A 
S, = | (Si X Sz) X Si Y = [A] [MCT] 
ai xe Ge Z Imes 
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<> > 
Where §, and Se are known unit vectors of stars in moon centered 
inertial space and MCI is the unknown coordinates of the spacecraft with 


respect to the moon centered inertial system. 


=> 
Sim xX 
> > ~ > 
Ss = (Sin X Sen)X Sin Y =  [B] [SM] 
Sim X Sem Z, } 


SM 


Where Sy and Sey are measured unit vectors of stars in stable 
member space of the stars, and SM is the assumed stable member 


coordinates referenced to the moon centered inertial coordinate system. 


[Sa] [Ss] 


[A] [MCT] [B] [SM] 


[SM] = [B]° [A] [Mc] 


[° REFSMMAT] = [B]° [A] 
Therefore: MCI = [{REFSMMAT]' [SM] 


Upon comparison with the gravity vector and an estimate of the lunar 


radius at the point, a position on the lunar surface can then be estimated. 


Considering only the precision of the counter reader (with interpolation) 
we See that the positional accuracy may be as poor as 152 meters 
or 498 feet. If we considered the errors in the determination of the unit 
Pravity vector, and the errors in the star tables, we would obviously find 
that this system is not adequate for determining accurate positions on the 
lunar surface. Considering a minimum standard error of 20",8, we will 


eliminate this system from consideration for geodetic purposes. 
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2.4.3 Method 3: Measurement of Relative Azimuth Angles 


Method three involves the simultaneous determination of longitude, 
latitude, and azimuth based on the motion of stars in azimuth. An 
ordinary theodolite can, of course be used, but it again suffers the 
drawbacks mentioned in Method 1. 

This method though, lends itself to a concept which can allow 
automatic determination of astronomic position and azimuth. An 
"Electronic Zenith Camera" (EZC) is currently under development by 
the Control Data Corporation of Minneapolis, Minnesota. It is described 
by Carroll as a portable sensor-computer combination which performs 
both detection of stars and calculation of astronomic position [1] [2]. 
Stars are detected by their transits across a fan of radial slits placed 
at the focal surface of a vertically oriented optical system. 

The principal advantages of this system is that once set up, the 
system can determine its astronomical position on the lunar surface 
without further attention from the astronauts. The system itself can 
be easily set up on the lunar surface. This set up includes approximate 
north orientation, and approximate levelling. If necessary, the system 
could be left to operate after the astronauts have left the lunar surface. 
Radio transmission of current updated positions could be relayed to 
earth if necessary. 

It is considered that astronomical longitude and latitude can be 
determined to within 0.3 (10) seconds and azimuth to within 3 (10) 
seconds (not considering star catalogue errors) which is certainly 
acceptable accuracy at this time. No instrument of this type has 
currently been manufactured, but all advance planning has been completed 
and the necessary technical skill is available to produce this instrument. 

The EZC can be broken down into the system block diagram as seen 


_— 


in Figure 5. 
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Figure 5. - EZC Block Diagram 


We will examine each component starting with the sensor 
(see Figure 6). The lens system is an f/3.3 system having a 3 
centimeter aperture and a 10 centimeter focal length. The lens system has 
a 16 degree total field of view (eight degrees maximum off-axis angle), and 
can form an image of a point source restricted to a tangential dimension 
no greater than 10 seconds of are but not restricted by the radial dimension. 
It will be necessary to shield the lens from such bright bodies as the sun 
and earth. 

The reticle has 100 radial slits (at the focal point of the lens) about 


1.4 centimeters long and 4.8 microns wide. 
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Figure 6. - EZC Sensor Unit 


21 





a 





Directly behind the slit reticle is placed an end-window photomultiplier 
tube which collects all light passing through the slits. The output of this 
tube is amplified, filtered, and processed so that only the maximum 
signal (peak pulse) will be recorded and timed. ‘The time input can be 
any commercially available temperature compensated crystal oscillator. 

Looking at the control unit, (see Figure 7) we find that the heart of 
the system is a miniature general purpose digital computer from which 
the sensor is controlled, and from which the operator receives status and 
calculated information. Also included are a precision oscillator (crystal 
oscillator), a power source, and the necessary interface electronics. The 
level reference is a mercury pool autocollimator. Miniature electronic 
autocollimators can reflect a beam off the surface of a eomnletely enclosed 
mercury pool (dust free surface), the surface of which is covered by a 
silicone fluid which damps out ripples. If the unit is to function during 
the lunar night, it will be necessary to heat the mercury pool, due to the 
lunar temperature being below the freezing point of mercury. 

The DC reversal motor is under command of the control unit. The 
motor drives the instrument about the vertical axis against one of two 
stops placed 180° apart. The reversal bearing is composed of dissimilar 
metals (brass and bronze) thereby eliminating the possibility of micro- 
scopic welds, which would inhibit accuracy. 

Prior to operation, the astronaut must level the instrument to within 
30 seconds. Within this range the automatic level reference supplies 
information to the computer so as to measure the angle between the 
instrument vertical axis: and the true local vertical. The astronaut must 
also make an initial north (south) - alignment of the instrument. This 
can be accomplished by using a gunsight mechanism on the equipment and 
sighting with the aid of near polar stars such as previously described. 
This will yield an initial approximation of north, with the exact azimuth 


to be solved for. © 
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We arc interested in detecting stars from +1. to +5.3 magnitude. 
It has been determined that the bcst results are obtained with a star 
transiting at least one slit each minute. In order to achieve this rate 
on the lunar surface, we must either wait approximately 27 times longer 
for a solution than we would on earth (approximately 27 hours) or have an 
internal spinning reticle which would reduce the solution time but would 
introduce additional errors into the system (approximately I-are second). 

Star identification is based on an approximate knowledge of the 
observer's position, and north alignment. The computer determines the 
expected pattern of transit times from the stored star catalog. This 
computed pattern and the measured time pattern are then superimposed, 
and these target pairs which correspond closely and uniquely in magnitude 
constitute tentative identifications. These are then used to adjust the 
assumed position and azimuth and thereby to adjust slightly the computer 
time pattern so as to obtain a better fit with the measured pattern. From 
here the final solution sequence preceeds by solving a linearized version 


of the following equation derived from equations 2.2: 


7) 
cot( Ao + B; ) sin(LHA, -A -Q@;,) -sin@cos (LHA, -r - a, 
+ tan dy; cos® = 0 


where 

Ao - The azimuth of the Oth slit (dependent on north alignment) 

B; - The fixed angle between the Jth slit and the Oth 
slit 

LHAG - Apparent lunar reference meridian hour angle of the node 
(LADT) 

r - Longitude 

Of - Lunar right ascension of the ith star 

D - Latitude 

dy - Lunar declination of the i th star 
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This linearized version is then 
0 = sin A; cot Z d® + cos ®( tan®- cosA; cot Z) dh - dA, 


where Z is the computed zenith distance 


Q 


h = LHAg -A-@, 


A; = Ao +B, 
and 


_ (sind sin9+ cos d_ cosh cos®) sin A, 
cou Z — 


~cosd sinh 


The measured quantities are azimuth and time. Latitude and longitude 


of the station are assumed. 


Since there are more measurements than unknowns, we may use a 
least square solution, with a Kalman filtering technique [10] based on the 
Newton-Raphson iteration formula, so that matrix inversions will not be 
required. This iterative procedure is based on the following equation: 


[ Ay Xy + Ly J 


XN4] = - Cy a 
ie NN ee eee ag 


where: 
d® | 9 | 
matrix Xy = sX: = |dh X, = 10] 
i | 
matrix ay = ia3 = (sin A; cot Zi - cos® (tan®-cosA; cot Z;)-1) 


2 
In = cot (A. +A;) sin(LHA, - XA -@; ) - sin® cos (LHA, 


~-rA-M@,) + tand; cos® 
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Cy an an Cr 
[1 + ay Cyan] 


and matrix Cy = Cy - 


I 


where matrix Ce 3C3 = ey) 


where the rows of A are the rows of ai to as. 


We can see that the incorporation of a new measurement does not 
require a new matrix inversion. Only one inversion is required and that 
is for Co which is only an inversion of a 3 x 3 matrix. After the first 
five marks no old observation information need be saved. ‘The values of 
the elements of matrix C are indicative of the quality of the results. 
Matrix C is the covariance matrix, but only under the assumption that 
there are no systematic errors in the residual vector L. There are 
systematic errors in the residual vector L, due to the timing error 
(7 libration ) and the errors in star right ascension and declinations. 

In order to cancel alignment errors in the optical-reticle system and 
normality errors between the instrument axis, the instrument is rotated 
ede and the astronomic position is again determined. The average of 
these positions is then used. The level status of the reversal bearing 


is also measured and it is applied as a final correction. 


2.4.4 Method 4: Measurement with Photography 


The determination of astronomical position by photography can be 
accomplished with a photo-theodolite or other metric camera using methods 
such as one or three. The drawbacks of a photo-theodolite are the same 
as for an ordinary theodolite and thus would tend to make it impractical 
for early missions. 

As described in the Geonautics study [7], a metric camera could be 


used as a zenith camera or otherwise fixed in elevation to photograph stars. 
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Various sophisticated systems such as gravity-operated automatic leveling 
devices with gimballed camera support, electronically operated automatic 
leveling devices, Syroscopic devices, etc., could be considered to level the 
camera. The simplest and most feasible device would seem to be the 
conventional spirit bubbles. These bubbles could be placed on a camera 
support such as a tripod, so that the camera, when placed on the support, 
could be theoretically leveled within the range of two orthogonally placed 


precise level vials inside the camera body. These precise levels can be 


photographed during each exposure, thus providing a permanent record of 
the level reference. Drawbacks in this system tend to lie in the fact that 
the bubbles are affected by differential heating from external light, so that 
they must be protected both from radiant heat and from light. The bubble, 
of the sensitivity required in the camera, would be termed a "'slow'' bubble 
in a conventional surveying instrument. It would be even slower on the 


lunar surface due to the moon's smaller gravitational force. 


This method provides a permanent and true record of field data as 
film can be removed from camera and returned to earth. Most advantageous 
is the fact that a great number of measurements can be taken in a minimum 
amount of time after the initial set up. There is a requirement for extra 
equipment such as a tripod, but it could be retractable and made of 
aluminum or some similar material Timing of each photo would also be 
required, but this could be done by an internal clock, which would furnish 


the time to be included in the edge data. 


2.5 Error Analysis 


Let us first consider the errors which contribute to the lunar apparent 
selenoequatorial coordinates. In the transformation of the mean geoequatorial 
coordinates to the mean ecliptic coordinates, the relatively small errors in 
a, 65, and € can be considered negligible. Therefore, 8 and A can be 


considered errorless. In the transformation of the mean ccliptic coordinates 
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to the mean sclenoecquatorial coordinates we consider that the error in 2 
is negligible, but there is definitely an error in inclination (I). Let us 
make a conservative estimate of 20" for this error. Differentiating 
equations (2.1) we see that the results are similar to that given by 


Mueller [16] . 
0 
ddz = cos@_ sinI d(A - &) + (sec dzsecBcosI - tan dz tan8) 
Q 
dp + sin@ , dI 


but d(A - 2) = 0 
ds = 0 


: (2 
therefore ddz = sin@® , dI and therefore the maximum value of 


Q Q 
wan sin@® , di = 20 0 sing ; 
82 2 ‘ . 
and dQ@, = (secd:cosI - secdstandz sinf) d (A - &) 
¢2 
-cos@, sec secd: sinId§ ~tandzcos@®, dI 


but again d(A -&) = 0 
dB = 0 


Q2 Q2 $2 
therefore d@, = -tand, cos@, dI = -20".0 tand: cos , 


Considering the transformation of mean into apparent selenoequatorial 


coordinates, we will use the following equations: 








c Q 0 Ai ar OA a1 0 Ai ax 
AQspp = d@, + — do + ST dil SL ddz 
OAi ay & OSBibi 0B: bi 0 Bi bi 
Say ip dk at iors dp + oes a dd, 
D (2.3) 
0 Bibi 9% <dCie1 8% dCi c SO Cr en 
5a iY d&@, + Baw d@, + 7, dI + AEE dda 
ODidi , S& SD, d, ~ OD: dy 
+ Sa d@,+t ST dI + 5 di. dda 
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and a similar expression for ddsapp, considering in both cases a 


negligible error in proper motion. 


In the following differential expressions we will assume: 


cos I 
dn 


0 
fe) 


1 


Limited to 45° 


1! 
1' 


The differential quantities will be assigned the following conservative 


values: 


ca = 
if = 
a 
dda = 
dT = 


Q 
d®@, = 


LO" 
Lo" 
20" 
20" 


15"! 
20" 


-000050 radians 


.000050 radians 
(estimate from [13] ) 


.000097 radians 


.000097 radians 
.000014 radians (estimate) 
~000097 radians 


Q2 
Ai a1 = Oe +0) (cosI - sinI sin@ , tandz ) 


3 Ai al 

— dog 
0 Al al 
a7 | : di 
OA al 
a 7 ddz 


= 0.0 


Ce 
(cosI - sinI sing 7 tan dz, )do 
Q 
(1 - .02676 sin@®, tandn)do = do 
LOO 


Q2 
(Pot +0 )(cosI sin@® , tand, - sinI)dI 
(-.0015059 ) (-. 97324) dI 
Or 0 


d . 8b 2 
(Pot +0)(sinI sin@® , sec dz) ddn 


(-.0015059) (. 05352) ddx 


29 





Q Q 
2A@ ag = (Plt +o) (~sinI cosas tands) da’ 
Oa G2 

7 = (-.0015059) (-.02676)d@, 
= o''.0 
Q 
Be ba = pcos@ , tandn 
3 Q | 
et ap = cos@ , tand;, dp . 
= 10.0 cosa tands — 
0 Bi bi 
a = 0".0 
on | st 
Q 
oa ddz =Pcos@_ sec’ dz dds 
= .00058 dda 
=O'.0 
Q $2 G2 
Bt a | = pein. tandg.da 
3 aS > ee 
= -(.00029) d@ x 
= QO'.0 
. | 2 Q 
Ci ci = K sin(lap - & + 180°) cosI cos@ , secdz dl 
Q 
at =~K sin (lap - &+ 180°) sinI cos@ ; secds dl 
= (.000097 )(1)(.02676)(1) (1.4) dI 
000 
, Q2 
= dd; = K sin (Lap - & + 180°) cosI tand: secdz cos@ g dda 


(. 000097) (1) (.02676) (1) (1) (1.4) dI 
0''.0 


I 
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3 a : a 
OO cy d@®., = Ksin (Lap - 8+ 180) Col cine’ 7 sced, aa 


dar? 


(.00097) (1) (.02676) (1) 1.4) aa 
0''.0 


(3 
K cos(Iap - 2 + 180°) sing, secds 
( 


Di dy = 
3 
SP ar = on.0 
ol 
OTh di .o) . 2 
sy dda = K cos (Lap - % + 180°) sin@ , tands secds dda 
Q 
= (000097) (1) (1) (1) (1.4) dda 
— 0"'.0 _ 
Tae &% 
9 Didi gq % =  Kesin (Lap - 8+ 180°) cos@®@, secdn d@, 
oes = (.000097) (1) dae 


= 06 


We will now consider the standard error in lunar right ascension. We 
will use the quadratic form of equation 2.3, eliminating terms which are 


zero. We will assume that dak and 


0 Bibi 
0p ee 
are dependent terms while all other non-zero terms are independent. 
Therefore, 
% 5 | = ay — bi 
(d@sapp) = ao’, jar (S doy” do)” 
fe) G2 


Q2 
400 tan’ dncos’ a 1 + 100 


+100 tan°dscos°a® 
- 200 tan*da cos’ a@ 


and ie = /300tan*dscos @, + 100 
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G2 
ho Bs = +(P5t Ol yes il ficou 


? 
0A} a} ne 


(2 
sinI cos@ , do 
80 


i 


Q2 
.02676 cos@nd0 


Il 


Max. value = 0".5 (will be neglected) 


3A1 at 2 
di =r +o )cosI cos@, dI 
= - .0015059 dI 
= aU 
OA ai Lo 
“Se dd; = 0".0 
dArart , tone 
don = -(Pit +o) -sinIsina, do, 
saY Q 
= (-.0015059) (-.02676) (1) (d@ a) 
; Q2 
Bi bi = psing y 


Q 
=9 Bibi g p= - sin@ , dpe 
dp 
’ 
phat = 0.0 


5 t 
CFP. dds = 0".0 


tL @ 
SBP G oe ee COs og cle 
On o 
=) =. 000Z0 de 
= O'.0 
t . Oo ‘ Q2 E 
Ci C1 = -K sin (Lap - +180) cosI sin@ » sindn 


-K sin (Lap - %%+ 180°) sinI cos dn 


32 





BO, C1 ; re) : : Q2 a 
sa, = (+K sin(Iap - $2 + 180 ) sinI sing , sin da 
~K sin(Lap - 21 + 180°) cosI cosd:) dI 
= -,.00007 dI 
=a 
: 02 
BS. ad, = (-Ksin (Lap - & + 180°) cosI sin@ , cos da 
pest 
+ K sin (Lap - & + 180°) sinI sind) dds 
= -,00007 dda 
= 0.0 
3 Pag C2 
sale ao n = -Ksin(Iap - &+ 180) cosIcos@, sindn 
Qn 
= .00007 aa 
= 0'.0 
' | ee Q. 
Di di = - K cos (lap - $2 + 180) cos@ 5 sind, 
5 t 
BS ar = 0". 9 
0D, dj ° Q2 
ad dd: = - Kcos(lap - $%@ + 180 ) cos@ , cos dg 
pest 
= .00007 dda 
= 0,0 
3 re QQ 02 
a das =  Keos (Lap - 8+ 180°) sin@ ; sindn 
Qa 
= .00007 daa? 
= O".0 
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We will now consider the standard error in lunar declination. We will 
use the quadratic form of the similar declination equation 2.3, eliminating 
terms which are zero. We will assume dd, and 

3 Bibi dp 
0p 
are dependent terms. 


Therefore, 8 Bibi d 2 
3p P) 





Gleam = (by 2 | 


’e) t 
A = Boda 
2,0 29 
= 400 sin @, + 100 sin Q@, 
Q2 
-200 sin® Qs 


2 
300 sin @, 


Q 
and ddsapp = 7 23 sind Q 


Thus, we can sce that if stars with near zero declination are used, we 
can effectively limit the error in lunar right ascension. Let us examine 

. fo) ee Oba 
the circumstances if we use stars with +10 declination and +10 right 


ascension. The maximum error will then be: 


dicen =\Vesou @oelooj jes 200 = 105 
tap = 1 ee SS 


Thus, we can see that there can be an optimum time for which to obtain 
a position. If we consider limiting gy to values of 350° right to 010° 
and 170° right to 190°, we will immediately limit the error in declination 
to a maximum value of 3.1. See Appendix E for stars which fall into 


this area. 
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Next we will consider errors in time. We have defined 


LADT 


LNMD i eed st 


or 
RADT 


II 


RMDT + eq.N 
We can consider that RMDT is errorless so that the error then occurs 
only in eq.N which is dependent upon o and T. So that we may say 
dT” = do? +dr® + dodtT 


or dT = 10°.0 draconic time or 5''.0 of are. 

We will now consider the standard error of each method in deter- 
mination of astronomic position. 

For Method 1, we will first present the differential formula for 


latitude [16] : 
d® = -secAdZ - cos@tanA dh 


- sin®cosec Z sec d - tand cot Z dd 
cos A 


We desire the azimuth to be near 0 or 180 degrees so that 


d2 = -dZ-(sin®cosec Z sec d - tan d cot Z) dd 


and if Z equals approximately 45°, then for a station at or near the equator 


d will also equal about 45° so that 
d’ = -dZ - dd 


Therefore, for Method 1, these errors are independent, letting dZ = 7"'.0, 


and considering 


(da) (7.0) + (17'.3 sinter 


and d® 


(SO etOm el 


Thus, an error in h will have little effect on latitude determination. The 


error in latitude determination is basically related to the approximate right 
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ascension. The differential for azimuth [16] is 


dA = (sin ®cosec”Z - sind cosec Z cot Z) dh 


+ cos ®sin h cosec*Z dd + cot Z sin A d® 
If we consider the zenith distance near 90° we see 
dA = (sin®- sind) dh + cos @sin hdd 
If the station is near the equator, © will be close to 0. Therefore, 
dA = - sin ddh + sin hdd 


Finally, we will select a pair of stars close to the meridian, preferably 
one with 90° north declination and one 90° south. Therefore, dA = 0 
for stations near the equator when we combine the pair. If we use only 
a sextant, we will have no means of accurately turning the angle between 
a star and landmark so that no permanent azimuth could be computed. 


The equation for longitude [16] is 


Q2 
= Chepeeot \7 acaieel Dea a Ti) 


where ae - the apparent right ascension 
- the hour angle as determined 
AL - the local time on time keeper 

AT, - the chronometer correction 


Il 


G2 
Therefore dA Bleeec., se Gln = CMD = l/NGE 


and dh = -(cosec A sec@dZ + cot A sec@d® 


-(tan®cosec Z sec d - tand cot Z sec d) dd 
sin A 
The error in h may be eliminated by observing at an azimuth of 90° or 
270° and by observing east and west stars in pairs which are symmetrical 
to the meridian and at the same altitude. Considering tan® is zero, the 


star pairs will eliminate the declination errors and zenith error terms. 
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02 
So that BP Oh EE fell lays 


d AT, has little effect on the moon for any ordinary chronometer, so that 


Q 
d= d@..- dT 


2 Q 
Therefore Sohe = (leas iP + (dT) - dTd@ uo 


2 2, 82 
= 300 tan~* d;2 cos @ , + 100 + 25.0 - 
eR Te ES 
5. Vaa00 tan” dz cos? + 100 


Limiting d:; to 45°, we see 
cA = 8''.6 to 18.0 


For Method 2, we have already stated that the low precision of the 
AOT precludes further consideration. 

For Method 3, with impersonal operation and high equipment precision, 
we can consider that the only errors presented will be systematic in nature, 
and these specifically will be the systematie errors in star coordinates 
and time. Thus, we can consider that an astronomic position by the use 


of this method will yield errors near the equator (declination < 10°) 


d®= ddap, = 3.1 (for optimum time) 


di= dof - at 


Q 
dX = (d@.,) +@T) - dT d@ap 
Ae Ote 
and dA = sin A, cot Z d+ cos (tan® - cosA, cot Z)dA 


near the equator, we say 


dA = (sin A, d®- cos A, dA) cotZ 


The largest Z we can hope for is 8°, and assigning an azimuth of the j™ 


slit to be 45° ,» we find: 


3f 





dA” = 25 d®? + 25d (considering no dependence) 


Therefore, dA. = 50" 30 


For Method 4, the errors are basically the same as for Method 1, 
if a panoramic camera is used, and basicaNy the same as for Method 3, 
if the zenith camera is used (considering camera operating at optimum 


time). See Table 3 for comparison of results. 


Method d& dr dA 
not directly 
BN oe foe Cn Oe low U 
: Bo eat obtainable 
se z not directly 
Z = 200 > 20u.8 Pee 
3 2 2 be@Pomaiuimy brie jo a2 > 50''.0 
4* Ss ae ney, SH) ERR I SR ERE 


*Dependent upon horizon or zenith camera type, not including 
instrument or observation errors. 


Table 3. - Estimated Standard Errors (Station Near the Equator) 


2.6 Conclusions 


Due to the possibility of limiting the error in star coordinates, it 
would seem that an automatic instrument as described in Method 3, would 
be the one which would be able to take best advantage of this point. 
Missions, of course cannot be scheduled so that these star coordinates 
could be utilized. It would be only coincidence if the current short 
duration landings occured at the proper time. Thus, the advantages of 


an instrument which could determine position automatically after the 


Soe 





astronauts depart is apparent. During future missions, if the landings are 
near the lunar equator, it would seem that the Electronic Zenith Camera/ 
Zenith photography would take best advantage of the stars with near zero 
declination. If landing were far removed from the equator, optimum 
times would need to be recomputed due to the increased error in.lunar 
right ascension, which would result from the higher star declinations. 
The format of the current missions, therefore, dictates that we use an 
automatic instrument which could take advantage of stars which offer 
minimum coordinate errors. The ease of set up and minimum training 
required for the astronauts, combined with the inherent accuracy of the 
system, show that the Electronic Zenith Camera or similar instrument is 
the best possible choice of equipment for the current missions. An 
accurate azimuth can only be determined by the use of an horizon camera 
at this time. 

It must again be stated that the Electronic Zenith Camera has not yet 


been manufactured, and could not be included in immediate missions. 
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3. DEFLECTION OF THE VERTICAL 


Deflections of the vertical relating the normals of the astronomic 


and selenographic systems can be obtained by: 


(1) Direct astronomical measurement of latitude and longitude 
with subsequent comparison with the selenodetic coordinates 
extended from a datum by triangulation or traverse. The 
results are affected by accumulation of error in the 


triangulation or traverse. 


(2) Integrating gravity anomalies over large areas surrounding 


the point where the deflection is required. 


(3) Computing the deflections of the vertical due to topography, 
and then refining this system using some model of isostatic 


compensation. 


(4) Computing the deflections of the vertical from known spherical 


harmonics. 


The first method depends on an assigned deflection at a datum, thus 
is of no use in this paper. The second method requires gravity anomalies 
over large areas, and these of course are not available except by obtaining 
them through method four. But, method four is of little value in deflection 
computations as the harmonics reflect moon-wide rather than local 
distribution. The computation of deflections of the vertical requires more 
local data near the station. Thus, we will investigate method three, but 
will use method four as a check on the magnitude of the deflections. 

In addition, the effect of other bodies on the deflection of the vertical 


will be calculated. 


3.1 Computation of Topographic Deflections of the Vertical 


Following the method of Jordan-Eggert [11] adapted for the moon, we 


consider the influence of the visible masses in the neighborhood of a station. 
AO. 








For this, we imagine the terrain around the station divided into compart- 
ments by means of concentric circles and radial lines, so that we obtain 
a great number of vertical columns between a given reference surface 
and the moon's physical surface. 

If the components of this deflection in the meridian and the prime 


vertical are expressed by AE and An respectively, then we have 


Seo : 
Aé = Rens (sin@z - sin@,) log, ro/ Qi 


t 


= Qt 
An a a (cos®2 - cos@,)logere/ri 
yt Fes] 


@ 
} 


where The surface layer density 


The mean density of the moon 


@ 


Height above reference surface 


oy 
i 


a - The radius of the moon 
AQ, Ae - The compartment azimuths 


Yi,Y2 - The outer and inner radius of compartments 


The reference surface which will be used is defined as a moon centered 


sphere of radius 1,737.988 kilometers. 


We will use the following values for the constants as given in [6]. 


S = 2.670 gm/cm* (assumed)* 
©, = 3.343 gm/cm* 

a = 1,737,988 meters 

0 = 206,265" 


—-= «+ 


* This figure was obtained from averaging fisures in various sources 
for a surface layer of 10KM. An error ia this figure will produce 
a proportional error in the deflections. 
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and we obtain 


AE = 0.0226 H(sin@. - sin@1) log. re/r, 
An = 0.0226 H(cos®@2- cos®,)log.er,/7r, 


Note: H in meters 


These two values are to be computed for each terrain section around 
the point 0, and then the two components of the deflection of the vertical 
for the point 0 are equal to the algebraic sum of all values of A€ or of 


An as: 


B n 
Cts ee n= & An, 
i=1 i=1 

We have so far assumed that the whole of the attracting stratum is 
in the horizon of the station, which for the moon would be grossly in 
error when only a short distance from the station. Therefore, it is 
necessary in order to secure the required degree of accuracy to take 
into account the curvature of the moon's surface. The effect of this 
curvature is to throw the compartment surface below the horizon of the 
station. Therefore, when significant, we will use the equation developed 


by A. R. Clarke, given in [11]. 


tan 1/4 Os 


AE=0". in@s - sin@ 
€ = 0'.0226 H (sin@2 - sind, ) log. tan 1/4 ©, 


+ cos 1/2 Os = cos 1/290, 


tan 1/4 O23 


An= - 0".02 Q2e- a 
n Zod (COsess.-) COs =) loa. any or 


+ cos 1/2 Os - cos 1/20, 
in which ©@=2 is re reduced to an are 
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in km 
re) - To in P ; 
7 1,737. 988 radians 


and ©, is ri reduced to an are 


S) = Yr; (in kn) 
i = 


1, 737.988 sale aT 


Any values for either the constant re / ri or for the constant difference 
(sin@2z - sin®,) may be adopted arbitrarily. Advantage is taken of this 
fact to adopt such values as would make the computation as simple and 
rapid as possible. The values re/ri = 1.19362053 and (sin®@g - sin@1) = 


0.25 were adopted. With these values, the deflection equations will be: 


Ag 


.001''h (in meters ) 


An 


.001" h (in meters) (the reference line is shifted 90°) 


As a result then, of this particular selection of arbitrary constants 
defining the limits of the compartments, the deflection produced at the 
station by the material lying above the reference surface in any compart- 
ment is expressed as one second of arc for each kilometer the mean 
elevation is above the reference surface. 

The arbitrary selection of the adopted values of (sin@2 - sin@, ) 
and re, ri was guided by three considerations. First, it was important 
to be as rapid as possible in determining the contribution of each section. 
The compartments must be small enough to bring the accuracy of the 
method within required limits and yet large enough to avoid excessive 
amounts of detail in the computation. The compartments should be 
compact areas, not long and narrow, in order to facilitate the estimation 
of the mean elevation within each compartment. 

The radii of the adopted circles separating the compartments are 
shown in Table 4. Each group of sixteen compartments (four in each 
quadrant) forms a ring. The rings, for convenience of designation, 
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were assigned serial numbers commencing with the inner ring, (see 
Figure 8). 

It was determined that the curvature of the moon had to be taken 
into account at ring 46. The outer radii of all succeeding rings were 
computed iteratively on the IBM 360-75 computer, (see Appendix A). 

It was found by Hayford [8], that on the earth, the computed 
topographic deflections are much larger than the actual deflections of the 
vertical, in that an influence is in operation which produces an incomplete 
counterbalancing of the deflections produced by topography, leaving much 
smaller deflections in the same direction. The deflections of the vertical 
are due to irregularities in the distribution of the masses composing the 
earth. These irregularities may occur either as a result of irregularities 
in the surface of the earth (topography), or as a result of irregularities in 
the distribution of the densities beneath the surface. 

The distribution of the densities below the surface of the earth is 
invisible and largely unknown, but the theory of Isostasy postulates pre- 
cisely such a relation between sub-surface densities and surface elevations. 

If the earth were composed of homogeneous material, its figure of 
equilibrium, under the influence of gravity and its own rotation, would be 
an ellipsoid of revolution. 

The earth though is composed of heterogeneous material which varies 
considerably in density. If this heterogeneous material were so arranged 
that its density at any point simply depended upon the depth of that point 
below the surface, or, more accurately, if all the material lying at each 
equipotential surface (rotation considered) was of one density, a state of 
equilibrium would exist and there would be no tendency toward a rearrange- 
ment of masses. 

If the heterogeneous material composing the earth were not arranged 
in this manner at the outset, the stresses produced by gravity would tend 
to bring about such an arrangement, but because the material is not a 


perfect fluid, the rearrangement will be imperfect. In the partial 
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RING Outer Radius (KM) RING Outer Radius (KM) 


1 . 0080 36 3.9211 
2 0095 ot 4.6803 
3 0114 38 De 0960 
4 - 0136 39 6.6682 
5) 0162 40 1.9592 
6 0194 Al 9.5003 
vi ~0231 42 11.3398 
8 0276 43 13.5354 
9 . 0330 44 16.1562 
10 0393 45 19. 2844 
11 . 0470 46 23.0184 
12 0561 AQ 27.4754 
13 . 0669 48 32. 1996 
14 20799 49 39.1461 
15 0953 50 46.7267 
16 - 1138 od Doe llot 
17 . 1358 a2 66.5781 
18 ~ 1621 Do 79.4741 
19 - 1935 04 94.8707 
20 2310 DO 113. 2545 
21 2157 26 135. 2082 
22 3291 ov 161.4304 
23 - 3928 28 192.7602 
24 -4688 a9 230. 2080 
20 -0996 60 274.9951 
26 . 6680 61 328.6057 
| 1973 62 392.8567 
28 ~9517 63 469.9958 
29 1.1359 64 062. 8462 
30 1.3558 65 675.0283 
31 1.6184 66 811.3203 
o2 dee le 67 ito ae ait 
Sh) 2.3057 68 1185.5146 
04 2.7922 69 1448. 1099 
S35) 3.2850 : 70 1793.1006 


Table 4. - Ring Outer Radius 
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rearrangement some stresses will still remain as different positions of 
the same horizontal stratum may have somewhat different densities, and 
the actual surface of the earth will depart from the ellipsoid of revolution 
in the sense that above each region of deficient density there will be a 
bulge or bump, and above each region of excessive density there will be a 
hollow. The bumps on the earth would be mountains and the hollows 
would be ocean floors. The excess material represented by that portion 
of the mountains above the reference surface will be compensated for by 
a defect of density in the underlying material and vice-versa for the 
oceans. This, of course, can be explained by the famous Pratt-Hayford 
theory of Isostasy. 

If we now consider the moon, we see that the consequences of 
hydrostatic theory as stated by Kopal [15] demand that given a sufficiently 
long time the material which on a short time-scale may behave as a 
solid is bound to get crushed under its own weight to settle to a form of 
minimum potential energy, whieh is a sphere. The departures of the 
lunar globe from spherical form permissible for a figure of equilibrium 
in the prevailing field of force are indeed small. To this extent, the 
moon may seem to be in hydrostatic equilibrium. However, this agree- 
ment becomes noticeably worse when we compare the momenta of the 
lunar globe with their theoretical values. Finally, by considering the 
mechanical ellipticity, we see that its value is so much at variance with 
the requirements of hydrostatic equilibrium that the discrepancy must be 
considered as real. 

Therefore, in spite of the approximate prevalence of hydrostatic 
equilibrium in the lunar interior, as attested by the nearly spherical form 
of its globe, the moon's motion around its center of gravity reveals the 
presence of small but unmistakable departures from this equilibrium. 

On the other hand, O'Keefe [19] was able to show by comparison between 
calculated Bouguer anomalies and the essentially free air anomalies of 


Muller and Sjogren, that the topography of the moon could be isostatically 
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compensated and that to a rough first approximation the moon is in 
isostatic equilibrium. He considers an actual crustal layer on the 
moon, for which he states the Pratt-Hayford theory of Isostasy could 
be applied. 

Under the assumption of Pratt's hypothesis of Isostasy, we will 
introduce a surface of compensation at the depth T below the reference 
surface. If we denote the changes of A€ and Ayn due to the changes in 
topography by MO&, and An, , then the sums A€ + AE, and An+ Ano 
are the terms which, when summed, give the total deflection of the 


vertical. From Jordan-Eggert [11]: 


| 2 
XE» = ao SO On py! (sinQ@s al sinQ@, )log. (. L T te ; 


47 ©x A r+ sae | 


8A@ Te" pga Ute 25 teem 
Eo = 47 @, A (cos@=o = cos @,)log, ee ee 
ey aN ie 
Aer = OH 


Therefore: 


2G) heme Ta. 
A&E, + AE = = Came cee sin®@,) [log. 7 

- log 72 a i - ra : ] 
ear eS 


VT? + pe? 
log, (8 
AE, + Aé 
let F = 35 = y -\4)4-—2-41 Vite +m 
lope = 
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Then for a definite value of T we can set up a table, from which for 
every value r, (because r2g = constant Xr,) the factor F can be 


taken. Then we will have 


AE + AE, = FAE 
An + An, = FAn 


O'Keefe states [19] that for 20% difference in lunar crust and lunar 
mantle density, a depth of compensation of 25 kilometers is in order, 
while for 10% difference, a depth of compensation of 50 kilometers is 
in order. Further extending these assumptions, calculations for depths 
of compensation of 25, 50, 75 and 100 kilometers has been made using 
the program and subprogram in Appendix A. 

Table 5 gives the various F factors for each ring owing to a 
specific assumed depth of compensation. 

Calculations have also been made for 125, 250, 500, 1,000 and 1,500 
kilometers depth of compensation primarily to obtain a full set of values. 
There is direct evidence of excess mass in the maria over what 

would be called for by Isostasy. The explanation of this excess mass 
which was discovered by Muller and Sjogren of JPL is subject to various 
interpretations. O'Keefe [19] states that this excess mass may be due 

to the maria being filled with some sort of material which arrived after 
the formation of the highland topography. H. C. Urey [21] states that 
this mascon effect could be produced by a flat circular slab of chrondritic 
material. For instance, Urey calculated that in Mare Imbrium, this 

disc would be some 4 kilometers thick and be 670 kilometers in diameter 
(mass = 5 X ay = grams). 

J. G. Stipe [20] states that these mascons could be caused by meteorites 
which have bombarded the moon and have been of necessary size and at 


necessary depth to create these mass anomalies. 
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Wise and Yates [22] state that mascons are produced by mantle plugs 
upwelling into giant impact bases punched through the lunar crust, followed 
by volcanic filling of the remainder of the crater above the plug. 

No matter which theory is accepted, the direct result will be a further 
deflection of the vertical. The area where the mascon is located, will 
deflect the vertical, proportionately larger than the height of this area 
would indicate. The exact numerical increase in the deflection would be 
dependent upon the mass distribution as offered in the theories above. 

We could consider these mass increases as additional point masses and 
compute their effect on the deflections accordingly, or we could select a 
station which would minimize their effect. 

Let us consider a central station which is quite distant from the major 
mass anomalies so that their effect will be insignificant (very low F factor) 
in the overall result. For instance, considering any station near the 
equator from 20° east longitude westward to 40° west longitude will limit 
the contribution of these mass anomalies on the deflection of the vertical 
as the station will be greater than 600 kilometers from the outer edge of 
these anomalies. We can see for a station greater than 600 kilometers 
from a large mass anomaly that the largest F factor is only about .015 
(down to 100 km compensation depth) and the contribution of the mass 
anomalies are then generally negligible on the overall result. For 
topographic deflections only, the area of these mass anomalies will 
generally fall into one compartment only, and will in turn affect the 
result very little. For this reason, I have selected the center of the 
crater Landsberg "A" as a starting point. Its coordinates are 328° 57'W 
and 0 13'N. This point was also selected due to availability of charts. 
Due to the nearly concentric area in terms of height of the center of 
Landsberg "A" and the use of a 1:250,000 chart, I have neglected rings 
1 to 35, which are approximately the inner 3 kilometers of rings. Rings 
06 to 50 were computed using the AMS series topographic lunar map sheet 


1 edition 1 of 3-64. This map is a modified stereographie projection 
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with a point of origin 4,410.10 kilometers from the point of tangency of 
the map plan and a selenographic sphere with a radius of 1,737.988 
kilometers. The contour interval is 250 meters based on a vertical 
datum centered at Mosting "A'™ at an elevation of 7,000 meters. The 
horizontal datum was also based on the center of Mosting ''A'' which was 
considered to be latitude 3° 10' 47" S and longitude 354° 50' 13", Rings 
51 to 63 were computed using an edition 1 - AMS sheet 2 topographic 
lunar map, scale 1:2,000,000 with the same datum, etc., as above. 
Rings 64 to 68 were computed using an edition 2 - AMS sheet 2 
topographic lunar map, scale 1:5,000,000 with the same datum, etc., 
as above. 

Rings beyond 68 were neglected due to their approach to the lunar 
limb; and their coinciding reduction in accuracy of height determination. 

Tables 6 and 7 give the heights of each compartment for computations 
in the prime vertical and in the meridian respectively. 

Table 8 shows the computed deflections for assumed depths of 


compensation. 
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Compensation Depth & n 


25 km + 2',934 + 10.983 

50 km + 3',332 + 20.709 
75 km + 2,234 O27 To 
100 km + 0.883 . + 33,123 
125 km - 0,547 + 37",447 
250 km - 6'',118 + 49.548 
500 kn - 16.792 +1'11'.984 
1,000 km - 30.184 +1' 38".947 
1,500 km - 38.529 a Agee 
= -1'04".450 | +2'17'', 480 


Table 8. - Deflections Considering Depths of Compensation 


We cannot consider that the value of © will remain constant down to 
the depths which we are considering. Thus, as depths increase, the 
average density for this surface layer will vary. We consider that the 
average density for these depths are as follows in Table 9 (interpolated 


from Kaula [12]). 
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KM Average Density Fil 


295 2.670 1 
20 3.000 1.1236 
79 3.140 1.1760 
100 3.210 1.2022 
125 3.200 1.2171 
150 3.263 1.2221 
250 3-336 1. 2494 
500 3.343 1.2521 
1,000 3.300 1.2360 
1,500 3.343 1.2521 


Table 9. - Average Surface Density for Selected Depth 


The F1 factor in Table 9 is a deflection multiplier, and its effect 


on deflections is given in Table 10. 
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Compensation Depth E n 


25 km + 2''. 934 + 10039383 
50 km coum Aurel io: ae P83) LO, 
To) kn + 2''.627 + 32'. 642 
100 km ‘ 1". 062 Too O20 
125 km = 0''. 666 + 45',447 
250 km S 7''. 644 Ply Oath O05 
500 km = ead" 5025 ted lueyans ibs ial 
1,000 km Seog a oUd ae aN a PAY 
1,500 km - 48" 242 +2'14''.509 
> - 1'04''. 450 +2'17".480 


Table 10. - Deflections Considering Changes in Density for 


Changes in Depth of Compensation 


The graph in Figure 10 shows the relationship of depth of compensation 
to deflections, considering no change in surface density, and considering a 


change in the surface density for changes in the depth of compensation. 
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Figure 9. - Deflections Considering Changes in Density for Changes in Depth of Compensation 





The maximum topographic deflections on the moon will most likely 
occur at a highland-maria dividing line. For instance, a north-south 
deflection was examined at the southern wall of the crater Pitatus at 
latitude 31°10' S and longitude 346° 10'. With a few assumptions (due 
to lack of large scale charts on the area) it was determined that this 
deflection was approximately 5 minutes. It is quite possible that other 
larger’ deflections exist, but probably will not greatly exceed this value. 
It is interesting to note that the deflections may be greatly minimized by 
choosing a station which is surrounded by nearly concentric rings of equal 
elevation. This will yield nearly zero deflections for the inner zones. 
This ideal station could be at the center of a concentric crater, at the 
top of an isolated rise on flat ground, or at a central point of a flat 
area. If, of course, a station were selected which did not satisfy this 


concentric condition, the inner rings must be considered. 


3.2 Computation of Spherical Harmonic Deflections of the Vertical 


In order to compare the values obtained by the method of isostatic- 
topographic deflections, we will use the formula as given by Dragg [3] to 


compute the deflections by the use of spherical harmonics. 





a0) 
| N Noa 
_ yD TRY) | = ing) - 
i EEO La (Aywcos MA + Byy sin MA) E | Pym (sind) 
' ON Eee), : 
sin® - N-1 Py +1, 2 (Sing) 
1 es esl 


M 
oe Yo ashy | 
MEG ee 5) Nou (Ayysin M XA - By; cos MA) Pyy (sin@) 





l= 


Dragg [3] neglected the 20, 21, and 40 spherical harmonic terms in 


his computation for the earth, but they have been included in this paper. 
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{ will use the spherical harmonics coefficients given by Blackshear [6], 


(see Appendix D). These are Cym and Sy coefficients which must be con- 
verted to Ay, and By above. . 
By definition we see that Cym = Jyw and Sy = Ky. 


For Heiskanen and Moritz [9], we see 
Aly = GM yr’ JNM 


Blym = -GMr Ky 


and that from Dragg [3] 


= 


( N - ji 
An = a Atm 
pv te 


Sl 
Bum ONS 


NM 
Nt+e2 


Thus we see that in terms of coefficients Ciyy and Sym 


eC ooUN 


gE = ae (CymcosMA + SywsinMA) [(N + 1) Pyw(sing) - 





sing - (N-M+1) Pysy m (sin) | 





7 = E uy M(CymsinMA - Sycos MA) Pym(sin®@) 
COS® 2 pO ) 


This equation was evaluated for the Landsberg "A" station using the 
IBM 360-75 computer, the program of which is contained in Appendix B. 
The results of this evaluation are 
€ = - 12'.184 
Qo 29138". 07 
These values fall between the isostatic-topographic results of 250 km 


depth of compensation, and the results of topographic compensation only, 
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when surface density changes are not considered. These values fall 
between the isostatic-topographic results of 250 kan and 1,500 km depths 


of compensation, when surface density changes are considered. 


3.3 Effect of Other Bodies on the Deflection of the Vertical 


It is obvious that the attractive effect of the sun and the earth will 
affect the deflection of the vertical on the lunar surface. The following 
developments refer directly to the earth but apply equally as well to the 
sun. Planets also have an effect but for all practical purposes this effect 
can be considered insignificant. 

Following the method of Jordan- Eggert [11] adapted for the moon, we 
see that 


3 Ne A’ sin Z cos Z 


Deflecti = 
ection M. e 
where M. - The mass of the earth 
M: - The mass of the moon 
A - The lunar radius 
Z - Zenith Distance of the earth 
r - Earth-moon distance 
Ae 
and where = §1.301 
Ma 
A” 1737.63) = 
4 _ (1787.68 _ g 
a ro (384, 400) eee 
Therefore: 
Deflection = 2.253 X 10°" sin Z cos Z 


I 


4". 646912274 sin Z cos Z 
Zes23 sin 2 Z 
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For the effect of the sun we will use the formula 


Son, A’ sin Z cos Z 





Deflection = Ma ¥ 
Where MoS 27,158,000 
Mn 
3 3 
A 1,737.63 fe 
ep en ot 1 BINS 3 10 
znd - (149,415, 834.4) 
Deflection = 0''.026 sin Z cos Z 


QueUorsraes 4 


which is negligible in comparison with other uncertainties introduced. 

We will now divide the deflection of the vertical caused by the earth's 
attraction into a northward and westward directed orponene: If A is 
the azimuth of the earth counted from the north, then we have for the two 


components of the deflection of the vertical 


fo — 21132 cin'2 2 cos A 
Nee". O2 Sin 2oZasin A 


These deflections are not constant for a station but vary in accordance 
with the change in earth-moon distance and the geometric librations in 
latitude and longitude. They also vary with the physical librations but 


these are insignificant. 


3.4 Conclusions 


Due to lack of knowledge of the lunar interior and its isostatic 
properties, we cannot, with any degree of confidence, state the magnitude 
of the deflections of the vertical at a given point. 

For the computation of the topographic deflections, we used lunar 


charts with a lunar centered spherical surface as the equipotential 
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reference surface, and from these the compartment relative heights were 
determined. If we consider the effect of an offset sphere on this computation, 
we can see immediately that the relative heights would change in accordance 
with the translation of the sphere's center. For instance, if the center 

of the reference sphere were offset toward the earth, the relative heights 
of the eastern compartments calculated in this paper would be consider- 
ably less, and the 7 deflections, considerably smaller. Thus, we will 
need better height determination on the moon, before the topographic 
deflections can be considered to be sufficiently accurate. Due to the 
relative ease of computation of the spherical harmonic deflections, it may 
be best to use them at present instead of the topographic deflections, 
especially in areas which have local concentric heights. The computed 
deflections could be, of course, corrected for the effect of the earth. 

If in the future, Method (1) is employed, and relative deflections 
between stations were obtained, we would be in a much better position to 
predict the absolute deflections at these same stations using the topo- 
graphic method. 

The employment of Method (2) would, if extensive enough, yield the 
absolute deflections of each station. It is considered that in this method, 
lies the best possible results. But, alternately, the most time and effort 
would be required over the other methods. 

Further work can be accomplished using Method (3). For instance, 
Karl Jung in Report No. 41 of The Ohio State University, Institute of 
Geodesy, Photogrammetry and Cartography, has given the formulas to 
obtain approximate values of the influence of lunar craters on the 
deflection. These computations could become particularly useful for 
stations offset from the crater centers, as they would yield a rapid 


evaluation of the deflection attributed to the inner rings of the template. 
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F FACTOR COMPUTATION PROGRAM AND ASSOCIATED 


MOON CURVATURE SUB- PROGRAM 


The main program computes the "F" factor for a given ring, given 


the depth of compensation. The associated sub-program computes the 


change in outer ring radius due to moon curvature. The main program list 


of variables is: 


IR 
J 
Al 
A2 
T 
FA 
FB 
FC 
FD 
FE 
FF 
FG 
Mal 
F 
The sub= 


Indexing variable for ring number 

Indexing variable and integer compensation depth 
Inner radius of ring 

Outer radius of ring 

Real number compensation depth 

Compensation depth squared plus outer radius squared 
Compensation depth squared plus inner radius squared 
Outer radius plus the square root of FA 

Inner radius plus the square root of FB 

FC divided by FD 

The natural log of FE 

The outer radius divided by the inner radius 

The natural log of FG 

The F factor equal to 1 - FF/FH 


program Mooncu uses an iterative procedure employing the 


secant method to compute the outer ring radius taking moon curvature into 


account. 


The sub-program list of variables is: 


THETAI 


THETA 2 


THETB2 


The inner radius divided by the lunar radius 
expressed in radians 

The outer radius uncorrected for lunar curvature 
expressed in radians 

An approximate outer radius corrected for 


curvature expressed in radians 
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Tl 
T2 
T2B 
T21B 
Alga! 
T11 
T3 
T3B 
SOC 


ASOC 


THETC2 


CSOC 
BSOC 


THETA1/4 

THETA 2/4 

THETB2/4 

Tangent of T2B 

Tangent of T2 

Tangent of T1 

ate all /-aleslal 

Wealiey ae 

The natural log of T3 minus the cosine of % 
THETA1 plus the cosine of $ THETA 2 minus 
the natural log of the outer radius divided by 
the inner radius which is given by definition. 
This is considered the first approximation. 
The natural log of T3B minus the cosine of $ 
THETAI1 plus the cosine of $THETB2 minus the 
natural log of the outer radius divided by the 
inner radius which is given by definition. 
This is considered the second approximation. 
By the secant method where: 


Pe ale.) nace de.CED, 
ti f(X1) ~ £(Xs-2) 


we see 


_ THETA2 (ASOC) - THETB2 (SOC) 
eae ASOC - SOC 


The difference between X,,, and X, 


Outer radius in kilometers 


fl 
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APPENDIX B 


LEGENDRE POLYNOMIAL COMPUTATION PROGRAM AND 


SPHERICAL HARMONIC DEFLECTION PROGRAM 


The program computes the Legendre polynomials and the deflections 


of the vertical in both the meridian and the prime vertical using the 13th 


degree polynomial coefficients as given by Blackshear [6] for the position 


of the center of Landsberg "A". The Legendre polynomials were 


computed to the fourteenth degree as required and the subscript 15 


represents the Oth element [P(15,15) = P(0,0)]. The program list of 


variables is: 


P(15, 15) 
C(14,14) 


(14, 14) 


ALAT 


ALONG 


I,J,KA,M 
MA,MB,N 
PSI 


APSI 
PNU 


APNU 


An array which contains the Legendre polynomials 
An array which contains the Blackshear C 
coefficients 

An array which contains the Blackshear § 
coefficients 

Latitude of computation point 

Longitude of computation point 

Sine of latitude 

Cosine of latitude 

Indexing variables 

Indexing variables 

Evaluation of a single term of the deflection 
in the meridian 

Sum of PSI Terms 

Evaluation of a single term of the deflection in 
the prime vertical 


Sum of PNU terms 


3 





PPSI Multiplication of APSI by constants for final 
value of deflection in meridian 
PPNU Multiplication of APNU by constants for final 


value of deflection in prime vertical 
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APPENDIX C 


PHYSICAL LIBRATION EQUATIONS 


As stated by Kolaczek [14], physical libration is the physical 
deviation of the moon's motion from that described by Cassini's Laws. 
It can be treated as perturbations on the elements of the lunar orbit 
and rotation. Thus, the three components are called libration in 
longitude 7, in node o, and in inclination Pp. The libration in 
longitude is a perturbation on the celestial longitude of the moon (can 
be considered as a correction to lunar time), libration in node on the 
celestial longitude of the node, and libration in inclination on the 
inclination I of the moon's equator. These perturbations are small, 


never exceeding 2' of arc. They can be expressed as follows: 


T = Daz sin (ng +mg' + pF +qD) 
Io = If +Xb: sin(ng +mg' + pF +qD) 
p = De, cos (ng + mg’ + pF +qD) 


Where g and g' are the mean anomaly of the moon and sun respec- 
tively (dependent on time). 
F = gtw 

D 


g-g'+w- wi 


Where w and w' are the arguments of the perigrees of the moon and 
sun (dependent on time), n, m, p, and q and the coefficients a, b, and 


c are given in Table 11. 


76 





Symbols Series Argument Multipliers Coefficients in 
n m p q Arc Seconds 


0 
0 
I 
Hl 
1 
1 
2 
2 
2 





Table 11. - Physical Libration Argument 


Multipliers and Coefficients 
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APPENDIX D 


BLACKSHEAR SPHERICAL HARMONIC COEFFICIENTS 


Table 12 below gives the C and S coefficients as determined by 
Blackshear in 1969. The coefficients are to the 13th order and degree. 


These coefficients are connected with satellite dynamics, with the 


potential V written in the form [9]: 
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APPENDIX E 


OPTIMUM POSITIONS FOR MINIMUM 
STAR COORDINATE ERRORS 


Stars which yield the minimum error in transformation were found 
to possess the coordinates of: 


0° 


Lunar Declination 


O 
0° or 180 


Lunar Right Ascension 


Thus this point will fall on the ecliptic at the position of the nodes. 
Figures 10 and 11 (from [16]) give the positions of these stars in earth 
coordinates of right ascension and declination for the time January 0.5 
1971, January 0.5 1972, January 0.5 1973 and January 0.5 1974. The 
motion of these points is in direct relationship with the motion of the 


lunar nodes. 
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